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Abstract. We show that the regularity of the gravitational metric tensor in spher- 
ically symmetric spacetimes cannot be lifted from C'^'^ to C^'^ within the class of 
C^'^ coordinate transformations in a neighborhood of a point of shock wave inter- 
action in General Relativity, without forcing the determinant of the metric tensor 
to vanish at the point of interaction. This is in contrast to Israel's Theorem [5] 
which states that such coordinate transformations always exist in a neighborhood 
of a point on a smooth single shock surface. The results thus imply that points of 
shock wave interaction represent a new kind of singularity for perfect fluids evolving 
in spacetime, singularities that make perfectly good sense physically, that can form 
from the evolution of smooth initial data, but at which the spacetime is not locally 
Minkowskian under any coordinate transformation. In particular, at such singu- 
larities, delta function sources in the second derivatives of the gravitational metric 
tensor exist in all coordinate systems of the C^'^ atlas, but due to cancelation, the 
curvature tensor remains uniformly bounded. 



1. Introduction 

The guiding principle in Albert Einstein's pursuit of general relativity was the prin- 
ciple that spacetime should be locally inertial 1. That is, an observer in freefall through 
a gravitational field should observe all of the physics of special relativity, except for 
the second order acceleration effects due to spacetime curvature (gravity). But the 
assumption that spacetime is locally inertial is equivalent to assuming the gravita- 
tional metric tensor g has a certain level of smoothness around every point. That is, 
the assumption that spacetime is locally inertial at a spacetime point p assumes the 
gravitational metric tensor g is smooth enough so that one can pursue the construction 
of Riemann Normal Coordinates at p, coordinates in which g is exactly the Minkowski 
metric at p, such that all first order derivatives of g vanish at p, with all second order 
derivatives of g bounded in a neighborhood of p. The nonzero second derivatives are 
then a measure of spacetime curvature. However, the Einstein equations are a system 
of partial differential equations (PDE's) for the metric tensor g and the PDE's by 
themselves determine the smoothness of the gravitational metric tensor by the evolu- 
tion they impose. Thus the condition on spacetime that it be locally inertial at every 



The ideas and methods presented here are the creation of Moritz Reintjes. The detailed proofs 
are due to him and can be found in his doctoral thesis '8', which was supervised by Blake Temple. 
Both authors were partially supported by second author's NSF Grant, where the problem was first 
proposed. 

^Also referred to as locally Lorentzian or locally Minkowskian. 
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point cannot be assumed at the start, but must be determined by regularity theorems 
for the Einstein equations. 

This issue becomes all the more interesting when the sources of matter and energy 
are modeled by a perfect fluid, and the resulting Einstein-Euler equations form a sys- 
tem of PDE's for the metric tensor g coupled to the density, velocity and pressure 
of the fluid. It is well known that the evolution of a perfect fluid governed by the 
compressible Euler equations leads to shock wave formation from smooth initial data 
whenever the flow is sufficiently compressive!! At a shock wave, the fluid density, pres- 
sure and velocity are discontinuous, and when such discontinuities are assumed to be 
the sources of spacetime curvature, the Einstein equations imply that the curvature 
must also become discontinuous at shocks. But discontinuous curvature by itself is not 
inconsistent with the assumption that spacetime be locally inertial. For example, if 
the gravitational metric tensor were C^'^, (differentiable with Lipschitz continuous first 
derivatives, [11]), then second derivatives of the metric are at worst discontinuous, and 
the metric has enough smoothness for there to exist coordinate transformations which 
transform g to the Minkowski metric at p, with zero derivatives at p, and bounded 
second derivatives as well, [H]. Furthermore, Israel's theorem, [5], (see also [H]) as- 
serts that a spacetime metric need only be C^'^, i.e., Lipschitz continuous, across a 
smooth single shock surface, in order that there exist a C^'^ coordinate transforma- 
tion that lifts the regularity of the gravitational metric one order to C^'^, and this 
again is smooth enough to ensure the existence of locally inertial coordinate frames 
at each point. In fact, when discontinuities in the fluid are present, C^'^ coordinate 
transformations are the natural atlas of transformations that are capable of lifting the 
regularity of the metric one order, while still preserving the weak formulation of the 
Einstein equations, [TT]. It is common in GR to assume the gravitational metric tensor 
is at least C^'^. For example, the C^'^ regularity of the gravitational metric is assumed 
at the start in singularity theorems of Hawking and Penrose, [4]. However, in Standard 
Schwarzschild Coordinate^ (SSC) the gravitational metric will be no smoother than 
C^'^, if a discontinuous energy momentum tensor in the Einstein equations is present, 

c.f. m- 

In this paper we prove there do not exist C^'^ coordinate transformations that can 
lift the regularity of a gravitational metric tensor from C^'^ to C^'^ at a point of a shock 
wave interaction in a spherically symmetric spacetime in GR, without forcing the de- 
terminant of the metric tensor to vanish at the point of interaction. This is in contrast 
to Israel's Theorem [5] which states that such coordinate transformations always exist 
in a neighborhood of a point on a smooth single shock surface. It follows that solu- 
tions of the Einstein equations containing single smooth shock surfaces can solve the 
Einstein equations strongly, (in fact, pointwise almost everywhere in Gaussian normal 
coordinates), but this fails to be the case at points of shock wave interaction, where 
the Einstein equations can only hold weakly in the sense of the theory of distributions. 
The results thus imply that points of shock wave interaction represent a new kind of 
singularity in General Relativity that can form from the evolution of smooth initial 



Since the Einstein curvature tensor G satisfies the identity Div G = 0, the Einstein equations 
G — kT imply Div T — 0, and so the assumption of a perfect fluid stress tensor T automatically 
implies the coupling of the Einstein equations to the compressible Euler equations Div T — 0. 

"^It is well known that a general spherically symmetric metric of form ds^ — ~A{t, r)dt^+B{t, r)dr'^ + 
E{t,r)dtdr + C{t,rfdQ,'^ can be transformed to SSC in a neighborhood of a point where ^ 7^ 0, c.f. 

m 



REGULARITY SINGULARITIES IN GENERAL RELATIVITY 



3 



data, that correctly reflects the physics of the equations, but at which the spacetime is 
not locally Minkowskian under any C^'^ coordinate transformation. At such singulari- 
ties, delta function sources in the second derivatives of the gravitational metric tensor 
exist in all coordinate systems of the C^'^ atlas, but due to cancelation, the curvature 
tensor remains uniformly bounded. 

To state the main result precisely, we consider spherically symmetric spacetime 
metrics g^u which solve the Einstein equations in SSC, where the metric takes the 
form 

ds^ = g^udx^'dx" = -A{t, r)dt^ + B{t, r)dr^ + r^dQ^, (1.1) 

where either t or r can be taken to be timelike, and dO,"^ = d{P + s\T?{'d)dip^ is the 
line element on the unit 2-sphere, c.f. |12] . In Section 2 we make precise the definition 
of a point of regular of shock wave interaction in SSC. Essentially, this is a point 
in (t, r)-space where two distinct shock waves enter or leave the point p at distinct 
speeds, such that the metric is Lipschitz continuous, the Rankine Hugoniot (RH) jump 
conditions hold across the shocks [ID], and the SSC Einstein equations hold weakly in 
a neighborhood of p and strongly away from the shocks. The main result of the paper 
is the following theorem, (c.f. Definition 13.11 and Theorem 17.11 below) : 

Theorem 1.1. Assume p is a point of regular shock wave interaction in SSC. Then 
there does not exist a C^'^ regular coordinate transformatioi^, defined in a neighborhood 
ofp, such that the metric components are functions of the new coordinates and such 
that the metric has a nonzero determinant at p. 

The proof of Theorem II. II is constructive in the sense that we characterize the Jaco- 
bians of coordinate transformations that smooth the components of the gravitational 
metric in a deleted neighborhood of a point p of regular shock wave interaction, and 
then prove that any such Jacobian must have a vanishing determinant at p itself. 
Because the metric becomes singular at p whenever regularity is forced upon it, 
we refer to points of regular shock wave interaction as regularity singularities. The 
numerical explorations in [14J strongly indicate that such singularities can form out 
of smooth initial data within a finite time just as with fluids governed by the special 
relativistic Euler equations, but we know of no complete mathematical proof of this 
fact. 

Our assumptions in Theorem 11.11 apply to the upper half [t > 0) and the lower 
half {t < 0) of a shock wave interaction (at t = 0) separately, suitable for the initial 
value problem, and also general enough to include the case of two timelik^ interacting 
shock waves of opposite familie^ that cross at the point p, but also general enough to 
include the cases of two outgoing shock waves created by the focusing of compressive 

^The atlas of C^'^ coordinate transformations is generic for lifting metric regularity at shock waves 
in General Relativity, because coordinate transformations preserve the metric regularity, (c.f. 
Section[5]), while a C^'" atlas seems to be appropriate only for metric tensors in C"'*^. The C^'^ atlas 
is natural because it preserves the weak formalism that derives from the quasilinear structure of the 
Einstein-Euler equations, a property we expect cannot be met by any atlas less regularity than C^'^, 
(e.g., a C°'^ atlas with resulting discontinuous metric components). Given this, points of regular shock 
wave interaction in SSC represent regularity stngulanties in the sense that they are points where the 
gravitational metric is less regular than in any coordinate system that can be reached within the 
C^'^ atlas. 

^In fact the theorem applies to non-null surfaces that can be regularly parameterized by the SSC 
time or radial variable, c.f. Theorem l7.1l b elow . 

*^Two shock waves typically change their speeds discontinuously at the point of interaction, c.f. [10] 
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rarefaction waves, or two incoming shock waves of the same family that interact at 
p to create an outgoing shock wave of the same family and an outgoing rarefaction 
wave of the opposite family, c.f. [10]. In particular, our framework and our theorems 
are general enough to incorporate and apply to the shock wave interaction which was 
numerically simulated in |14] . 

Historically, the issue of the smoothness of the gravitational metric tensor across 
interfaces began with the C^'^ matching of the interior Schwarzschild solution to the 
vacuum across an interface, followed by the celebrated work of Oppenheimer and Sny- 
der [9] who gave the first dynamical model of gravitational collapse. In |llj . Smoller 
and Temple extended the Oppenheimer-Snyder model to a C°'^ shock wave interface 
that allowed for nonzero pressure. In his celebrated 1966 paper [5], Israel gave the 
definitive conditions for C^'^ matching of gravitational metrics at smooth interfaces, 
by showing that if the second fundamental form is continuous across a single smooth 
interface, then the RH jump conditions also hold, and Gaussian normal coordinates 
smooth the metric to C^'^ and thereby provides a locally inertial coordinate system. 
In [3] Groah and Temple addressed these issues rigorously in the first general existence 
theory for shock wave solutions of the Einstein-Euler equations. In coordinates where 
their analysis is feasible, SSC, the gravitational metric was Lipschitz continuous at 
shock waves, but could be no smoother, and it has remained an open problem as to 
whether the weak solutions constructed by Groah and Temple could be smoothed to 
C^'^ by coordinate transformation, like the single shock surfaces addressed by Israel. 
The results in this paper resolve this issue by proving definitively that the weak solu- 
tions constructed by Temple and Groah cannot be smoothed within the class of C^'^ 
coordinate transformations when they contain points of shock wave interaction. 

As a final comment, we note that although points of shock wave interaction are 
straightforward to construct for the relativistic compressible Euler equations in flat 
spacetime, to our knowledge there is no rigorous mathematical construction of an 
exact solution of the Einstein equations containing a point of shock wave interaction, 
where two shock waves cross in spacetime. But all the evidence indicates that points 
of shock wave interaction exist, have the structure we assume in SSC, and in fact 
cannot be avoided in solutions consisting of, say, an outgoing spherical shock wave 
(the blast wave of an an explosion) evolving inside an incoming spherical shock wave 
(the leading edge of an implosion). The existence theory of Temple and Groah [3j 
lends strong support to this claim by establishing existence of weak solutions of the 
Einstein-Euler equations in spherically symmetric spacetimes. The theory applies to 
arbitrary numbers of initial shock waves of arbitrary strength, existence is established 
beyond the point of shock wave interaction, and the regularity assumptions of our 
theorem are within the regularity class to which the Groah- Temple theory applies. 
All of this is substantiated by the recent work of Vogler and Temple which gives a 
numerical simulation of a class of solutions in which two shock waves emerge from a 
point of interaction, and the numerics demonstrate that the structure of the emerging 
shock waves meets all of the assumptions of our theorem at the point of interaction. 
Taken on whole, we interpret this as definitive physical proof that points of shock wave 
interaction exist in GR, and meet the regularity assumptions of our theorem. Given 
this, the conclusion of our theorem is that points must exist where the gravitational 
metric tensor cannot be smoothed from C^'^ to C^'^ by coordinate transformation. 

In Section [2] we begin with preliminaries. In section [3] we set out the framework of 
shock waves in GR, and define what we call a point of regular shock wave interaction 
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in SSC. In Section[3]we give a precise sense in which a function is said to be C*^'^ across 
a hypersurface, and we introduce a canonical form for such functions, which isolates 
the Lipschitz regularity from the regularity in a neighborhood of the hypersurface. 

In Section [5] we derive conditions on the Jacobians of general C^'^ coordinate trans- 
formations necessary and sufficient to lift the regularity of a metric tensor from C"'^ 
to at points on a shock surface. The result is a canonical form for the Jacobians of 
all coordinate transformations that can possibly lift the regularity of the gravitational 
metric tensor to C^. 

In Section [6l we give a new constructive proof of Israel's theorem for spherically 
symmetric spacetimes, by showing directly that the Jacobians expressed in our canon- 
ical form do indeed smooth the gravitational metric to C^'^ at points on a single 
shock surface. The essential difficulty is to prove that the freedom to add an arbi- 
trary C^-function to our canonical form, is sufficient to guarantee that we can meet 
the integrability condition on the Jacobian required to integrate it up to an actual 
coordinate system. The main point is that this is achievable within the required 
gauge freedom if and only if the RH jump conditions and the Einstein equations hold 
at the shock interface, |10] . 

The main step towards Theorem 11.11 is then achieved in Section [7] where we prove 
that at a point of regular shock wave interaction in SSC there exists no coordinate 
transformations of the (t,r)-plane that lift the metric regularity to C^. The essential 
point is that the gauge freedom in our canonical forms cannot satisfy the inte- 
grability condition on the Jacobians, without forcing the determinant of the Jacobian 
to vanish at the point of interaction. In section [8] we extend this result to the full 
spacetime atlas of coordinate transformations that allow changes of angular variables 
as well as thereby proving Theorem ll.il In the final section [H we we show that 

Theorem Theorem 11.11 implies the non-existence of locally inertial frames at points of 
regular shock wave interaction. Since we do not know how to make mathematical and 
physical sense of coordinate transformations less regular than C^'^ in general relativity, 
we conclude that points of shock wave interaction represent a new kind of regularity 
singularity in spacetime. 



2. Preliminaries 

In General Relativity, the gravitational field is described by a Lorentzian metric 
g of signature (—1,1,1,1) on a four dimensional spacetime manifold M. We call 
M a C'^-manifold if it is endowed with a C'^-atlas, a collection of four dimensional 
local diffeomorphisms that are regular from M to M^. A composition of two local 
diffeomorphisms x and y of the form xoy~^ is refered to as a coordinate transformation. 
In this paper we consider C^'^-manifolds, since this low level of regularity offers a 
generic framework to address shock wave solutions of the Einstein-Euler equations. 

We use standard index notation for tensors whereby indices determine the coordinate 
system, (e.g., Tlf denotes a (l,l)-tensor in coordinates and Tj denotes the same 
tensor in coordinates x^), and repeated up-down indices are assumed summed from 
to 3. Under coordinate transformation, tensors transform by contraction with the 
Jacobian 



(2.1) 
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Ji denotes the inverse Jacobian, and indices are raised and lowered with the metric 
and its inverse g"^^ , which transform as bilinear forms, 

g^.u = JlA9^3^ (2.2) 

c.f. |15j . We use the fact that a matrix of functions Jj^ is the Jacobian of a regular 
local coordinate transformation if and only if the curls vanish, i.e., 

= Jj;. and Det ( Jf) / 0, (2.3) 

where fj = denotes partial differentiation with respect to the coordinate and 
Det (^jf^ denotes the determinant of the Jacobian. 

The time evolution of a gravitational field in general relativity is governed by the 
Einstein equations [T] 

G'^ + Ag'^ = kT'^ , (2.4) 

a system of 10 second order partial differential equations that relate the metric tensor 
gij to the undifferentiated sources T^^ through the Einstein curvature tensor 

G'^ = R'^ - ^Rg'^ , (2.5) 

a tensor involving second derivatives of g. Here A is the cosmological constant, (our 
results apply to a vanishing as well as a non- vanishing cosmological constant), k = 
is the coupling constant which incorporates Newton's gravitational constant G and the 
speed of light c, and T^^ is the stress energy tensor. 

We assume throughout that T*-^ is the stress tensor for a perfect fluid, 

= {p + p)uiUj+pgij, (2.6) 

where p is the energy density, Ui the 4-velocity, and p the pressure. Conservation of 
energy and momentum enter the Einstein equations through, 

T'{j=0. (2.7) 

which reduces to the relativistic compressible Euler equations in flat spacetime, and 
follows from the divergence free property of the Einstein equations, G*-'.^- = a property 
built into G at the start as an identity following from the Bianchi identities of geometry. 
Here as usual, semicolon denotes covariant differentiation v-j = vj +r^^t;', where F^^- 
denote the Christoffel symbols or connection coefficients associated with metric g, c.f. 

m- 

System ()2.4p and (j2.7p forms the coupled Einstein-Euler equations, a system of 
second order differential equations for the unknown metric gij, and unknown fluid 
variables p, p and . For example, imposing an equation of state p = p{p) closes the 
system, yielding a set of fourteen differential equations in fourteen unknowns. In special 
relativity the spacetime metric is taken to be gij = rjij where r]ij = diag{— 1,1,1,1) 
is the Minkowski metric, in which case (j2.7|) reduces to the relativistic compressible 
Euler equations, a system of conservation laws in which it is well known that shock 
waves form out of smooth initial data whenever the flow is sufficiently compressive, 
(see |10] or [7]). Shock waves are discontinuous solutions that only solve the Euler 
equations weakly, in a distributional sense. Across a smooth shock surface S, the 
Rankine-Hugoniot jump conditions hold, 

[T^nu = 0, (2.8) 
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where [/] = fi — Ir denotes the jump from right to left (wrt r) in function / across S, 
and riy is the surface normal. In particular, for smooth shock surfaces, the jump con- 
ditions (j2.8p are equivalent to the shock wave solution satisfying the weak formulation 
of dlZD across S, c.f. [lOj. 

For many astrophysical processes it makes sense to assume the spacetime is spher- 
ically symmetric, by which we mean a Lorentz manifold admitting two Killing vector 
fields that give rise to coordinates in which the metric takes the simplified form 

ds'^ = -Adt^ + Bdr'^ + 2Edtdr + CdO^ , (2.9) 

where the coefficients A, B, C and E depend only on t and r, and dfi^ is the line 
element on the unit 2-sphere, |15j . In this case we often suppress the dependence 
on (p and 9, and refer to the spacetime parameterized by the variables t and r as the 
{t, r)-plane. 

In a spherically symmetric spacetime with ^ ^ 0, one can always transform to 
Standard Schwarzschild Coordinates (SSC), where the metric takes the form [15j, 

ds^ = -A{t, r)dt^ + B{t, r)dr'^ + r^dO^. (2.10) 

The Einstein equations for a metric in SSC are given by 

B^ + B^^ = kAB^vT^^ (2.11) 
r 

Bt = -nAB'^rT^^ (2.12) 

Ar-A^-^ = kAB^vT^^ (2.13) 
r 

Bu-Arr + <^ = -2KABr^T^'^ , (2.14) 

with 

BAtBt Bf Ar ABr Al ArBr 
(J) = — * L _| L _| L _| 1 L . 

2AB 2B r rB 2A 2B 
The first three Einstein equations in SSC play a crucial role in the method we develop 

in this paper, and it is straightforward to read off from the first three equations that 
the metric cannot be any smoother than Lipschitz continuous if the matter source is 
discontinuous, for example, T^^ € L°°, and in this paper we make the assumption 
throughout that the gravitational metric is Lipschitz continuous. This provides a 
consistent framework to address shock waves in GR, in agreement with the theory 
and examples of shock wave solutions to the coupled Einstein-Euler equations, see for 
instance [H] or [3]. 

Lipschitz continuity arises naturally in the general problem of matching two space- 
times across a hypersurface, as first considered by Israel in [5j. Israel proved the rather 
remarkable result that whenever a metric is Lipschitz continuous {C^'^) across a smooth 
single shock surface S, there always exists a coordinate transformation defined in a 
neighborhood of S, that smooths the components of the gravitational metric to C^'^. 
The precise result is that the gravitational metric is smoothed to C^'^ in Gaussian 
Normal Coordinates if and only if the second fundamental form of the metric is con- 
tinuous across the surface. The latter is an invariant condition meaningful for metrics 
Lipschitz continuous across a hypersurface, and is often referred to in the literature as 
the junction condition, c.f. [15]. In [TT], Smoller and Temple showed that in spherically 
symmetric spacetimes, the junction conditions hold across radial surfaces if and only 
if the single [T*-^] UiUj = 0, implied by (j2.8p . holds. Thus, for example, single radial 
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shock surfaces can be no smoother than Lipschitz continuous in SSC coordinates, but 
can be smoothed to C^'^ by coordinate transformation. However, it has remained an 
open problem whether or not such a theorem apphes to the more comphcated C^'^ SSC 
solutions proven to exist [3]. Our purpose here is to show that such solutions cannot 
be smoothed to in a neighborhood of a point of regular shock wave interaction, a 
notion we now make precise. 

3. A POINT OF REGULAR SHOCK WAVE INTERACTION IN SSC 

In this paper we restrict attention to radial shock waves, hypersurfaces S locally 
parameterized by 

j:{t,^,v) = {t,x{t),'d,^), (3.1) 

and across which T in (j2.6p is discontinuous. For such hypersurfaces in SSC, the 
angular variables play a passive role, and the essential issue regarding smoothing the 
metric components by C^'^ coordinate transformation lies in the atlas of coordinate 
transformations acting on the {t, r)-plane alone. (The main issue is to prove theorems 
16.11 and 17.11 for (t, r)-transformations. In section [8] we discuss the straightforward 
extension to the full atlas of transformations that include the angular variables.) Thus 
we introduce 7, the restriction of a shock surface S to the {t, r)-plane, 

7(t) = (t,x(t)), (3.2) 

with normal 1-form 

n„ = {x,-l). (3.3) 

For radial shock surfaces (|3.ip in SSC, the RH jump conditions (j2.8p take the simplified 
form 

[tOO]x = [tOI], (3.4) 
[riO]x = [Til]. (3.5) 
Now suppose two timelike shock surfaces Sj are parameterized in SSC by 

j:i{t,e,<j)) = {t,x,it),e,cj)) i = i,2. (3.6) 

Let ji{t) denote their corresponding restrictions to the (t,r)-plane, 

7,(t) = (t,Xi(t)), (3.7) 

with normal 1-forms 

{ni)„ = {xi,-l), (3.8) 

and use the notation that [f]i{t) = [/(7i(i))] denotes the jump in the quantity / across 
the surface 7i(t). 

For our theorem it suffices to restrict attention to the lower or upper part of a shock 
wave interaction that occurs at t = 0. That is, in either the lower or upper half plane 

Rl = {{t,r) :t<0}, 

or 

M.l = {{t,r):t>0}, 



^Note that if t is timelike, then all timelike shock surfaces in SSC can be so parameterized. Our 
subsequent methods apply to spacelike and timelike surfaces alike, (inside or outside a black hole, c.f. 
[12]), in the sense that t can be timelike or spacelike, but without loss of generality and for ease of 
notation, in the remainder of this paper we restrict to timelike surfaces parameterized as in (|3.ip . 
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respectively, whichever half plane contains two shock waves that intersect at p with 
distinct speeds. (We denote with the closure of M^.) Thus, without loss of gener- 
ality, let 7j(t) = (t,Xi{t)), {i = 1,2), be two shock curves in the lower (t,r)-plane that 
intersect at a point (0,ro), tq > 0, i.e. 

a;i(0) =ro = X2(0). (3.9) 

With this notation, we can now give a precise definition of what we call a point of 
regular shock wave interaction in SSC. By this we mean a point p where two distinct 
shock waves enter or leave the point p with distinct speeds. The structure makes pre- 
cise what one would generically expect, namely, that the metric is a smooth solution of 
the Einstein equations away from the interacting shock curves, the metric is Lipschitz 
continuous and the RH jump conditions hold across each shock curve, and derivatives 
are continuous up to the boundary on either side. In particular, the definition refiects 
the regularity of shock wave solutions of the coupled Einstein-Euler equations consis- 
tent with the theory in ^ and confirmed by the numerical simulation in [14J. Without 
loss of generality we assume a lower shock wave interaction in M^. 

Definition 3.1. Let tq > 0, and let g^y he an SSC metric in C°'^ (^A/'nM^^, where 

A/^ C is a neighborhood of a point p = (0,ro) of intersection of two timelike shock 

curves ^dt) = it^xAt)) G R^, t G (— e, 0). Assume the shock speeds Xj(0) = limxj(t) 

t— s-O 

exist and are distinct, and let J\f denote the neighborhood consisting of all points in 
Af n ]R?_ not in the closure of the two intersecting curves 7i(t). Then we say that p is 
a point of regular shock wave interaction in SSC if: 

(i) The pair {g,T) is a strong solution of the SSC Einstein equations (j2.1ip - (l2.14p 

in N, with T^'" G C°iAf) and g^^ G C^iM). 
(a) The limits of T and of metric derivatives g^u,a exist on both sides of each shock 

curve jiit) for all —e<t<0. 
(Hi) The jumps in the metric derivatives [g^u,cT]i(t) are function with respect to t 

for i = 1, 2 and for t G (— e, 0). 
(iv) The limits 

limb/.i.,a]i(i) = [g^u,a]i{0) 

t-¥0 

exist for i = 1,2. 

(v) The metric g is continuous across each shock curve ")i{t) separately, but no better 
than Lipschitz continuous in the sense that, for each i there exists n, v such that 

at each point on 7^, t G (— e, 0) and 

\\vci\g^y^a\i{niT ^ 0. 

(vi) The stress tensor T is bounded on Af (1 M?_ and satisfies the RH jump conditions 

[r^"Mni)^ = o 

at each point on 7i(t), i = 1, 2, f G (— e, 0), and the limits of these jumps exist up 
to p as t 0. 
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4. Functions C°'^ across a hypersurface 

In this section we give a precise definition which isolates the C*^'^ regularity of a 
function to a hypersurface, and use this to derive a canonical form for such functions. 

Definition 4.1. Let Tj be a smooth (timelike) hypersurface in some open set M C M'^. 
We call a function f "Lipschitz continuous across S ( or C^'^ across T,), if f € 
C^'^Af), f is smootlE in A/" \ S, and limits of derivatives of f exist and are smooth 
functions on each side of S separately. We call a metric g^p Lipschitz continuous 
across S in coordinates if all metric components are across n. 

The main point of the above definition is that we assume smoothness of /, (or g^;/), 
away and tangential to the hypersurface S. Note that the continuity of / across S 
implies the continuity of all derivatives of / tangent to S, i.e., 

1/,.]^''^ = 0, (4.1) 

for all tangent to S. Moreover, Definition 14.11 allows for the normal derivative of / 
to be discontinuous, that is, 

\IM + 0, (4.2) 
where n'^ is normal to E with respect to some (Lorentz-) metric g^^ defined on M . 

We can now clarify the connections between the Einstein equations and the RH jump 
conditions (|3.4p . (j3.5p for SSC metrics only C*^'^ across a hypersurface. So consider 
a spherically symmetric spacetime metric (jl.ip given in SSC, assume that the first 
three Einstein equations (|2.1ip - (j2.13p hold, and assume that the stress tensor T is 
discontinuous across a smooth radial shock surface described in the (i, r)-plane by 
7(t) as in (|3.ip - (|3.3p . To this end, condition (|4.ip across 7 applied to each metric 
component g^.^ in SSC ()2.10p reads 

\Bt\ = -x\Br\ (4.3) 

\At\ = -x[Ar\. (4.4) 

On the other hand, the first three Einstein equations in SSC ()2.1ip - ()2.13p imply 

[B,\ = KAB'^r[T'^\ (4.5) 

[Bt] = -K^i?V[TO^], (4.6) 

[Ar] = KAB'^r[T^\ (4.7) 

Now, using the jumps in Einstein equations (|4.5p - (j4.7p . we find that (j4.3p is equivalent 
to the first RH jump condition (j3.4p l^ while the second condition (14. 4p is independent 
of equations (I4.5p - (l4.7p . because At does not appear in the first order SSC equations 
(I2.1ip -( i2.13p . The result, then, is that in addition to the assumption that the metric 
be C'^'^ across the shock surface in SSC, the RH conditions (13. 4j) and (13. 5j) together 
with the Einstein equations (I4.5j) - (14.7p . yield only one additional condition over and 
above (|4.3j) and (14. 4p . namely, 

[A,] = -x[Bt] . (4.8) 
The RH jump conditions together with the Einstein equations will enter our method 
in Section [5113 only through the three equations ()4.8p , ()4.3p and ()4.4p . 



For us, "smooth" means enough continuous derivatives so that smoothness is not an issue. Thus 
here, / £ C^(A/'\ E) suffices. 

'^This observation is consistent with Lemma 9, page 286, of [11], where only one jump condition 
need be imposed to meet the full RH relations. 
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The following lemma provides a canonical form for any function / that is Lipschitz 
continuous across a single shock curve 7 in the {t, r)-plane, under the assumption that 
the vector n^, normal to 7, is obtained by raising the index in (j3.3p with respect to a 
Lorentzian metric g that is C^'^ across 7. (Note that by Definition 14.11 varies 
in directions tangent to 7, and we suppress the angular coordinates.) 

Lemma 4.2. Suppose f is C^'^ across a smooth curve j{t) = (t,x{t)) in the sense of 
Definition \4-l\ t S (— e,e), in an open subset Af o/M^. Then there exists a function 
$ e C^(7V) such that 

f(t,r) = ^ip{t)\x{t)-r\+^t,r), (4.9) 

where 

^(i) = %l^ G C\-e,e), (4.10) 

and n^{t) = (i;(t),— 1) is a 1-form normal to the tangent vector v^{t) = j^{t). In 
particular, it suffices that indices are raised and lowered by a Lorentzian metric g^^y 
which is C^'^ across 7. 

In words, the canonical form (|4.9p separates off the C'''^ kink of / across 7 from 
its more regular behavior away from 7: The kink is incorporated into \x{t) — r[, 
(/9 gives the smoothly varying strength of the jump, and $ encodes the remaining 
behavior of /. I.e., 99 gives the strength of the jump because upon taking the jump 
in the normal derivative of / across 7, the dependence on cancels out. Note finally 
that the regularity assumption on the metric across 7 is required for ip{t) to be well 
defined in (j4.10p . and also to get the regularity of n^Ufj tangent to 7. In Section 
[6] below we prove Israel's Theorem for a single shock surface by constructing a C^'^ 
coordinate transformation using lTOl of Lemma l4.2l as a canonical form for the Jacobian 
derivatives of the transformation. 

In Section [7] we need a canonical form analogous to (14. 9p for two shock curves, but 
such that it allows for the Jacobian to be in the weaker regularity class C^'^ away from 
the shock curves. To this end, suppose two timelike shock surfaces described in the 
(t, r)-plane by, 7j(t), such that (I3.6p - (l3.9p applies. To cover the generic case of shock 
wave interaction, we assume each 7i(t) is smooth away from t = with continuous 
tangent vectors up to t = 0, and it suffices to restrict to lower shock wave interactions 
in 

Corollary 4.3. Let ji{t) = {t,Xi{t)) be two smooth curves defined on I = (— e,0), 

some e > 0, such that the limits lim 7i(t) = (0,ro) and Xi(0) = lim Xi{t) both exist 

t^o- t— >-0- 

for i = 1,2. Let f be a function in C^'^{M PI M?_) for M a neighborhood of (0, rg) in 

M^, so that f meets condition on each 7,. Then there exists a C^'^ function <I> 

defined on Af H M?_, such that 

= = i = l,2, (4.11) 

and 



2 

j=l,2 



for all [t, r) in J\f H M?_, where 



^iit) = i^^^ G C'>^\I). (4.13) 
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In particular, ipi has discontinuous derivatives wherever /o^j does, and again it suffices 
that indices are raised and lowered by a Lorentzian metric g^y which is C^'"^ across each 

5. A Necessary and Sufficient Condition for Smoothing Metrics 

In this section we derive a necessary and sufficient pointwise condition on the Jaco- 
bians of a coordinate transformation that it lift the regularity of a C^'^ metric tensor 
to C^'^ in a neighborhood of a point on a single shock surface S. In the next section 
we use this condition to prove that such transformations exist in a neighborhood of a 
point on a single shock surface, and in the section following that we use this pointwise 
condition on each of two intersecting shock surfaces to prove that no such coordinate 
transformation exists in a neighborhood of a point of shock wave interaction. 

We begin with the transformation law 

9aP = Ja9fMuJp, (5.1) 

for the metric components at a point on a hypersurface S for a general C^'^ coordinate 
transformation x", where, as customary, the indices indicate the coordinate 

system. Let Ja denote the Jacobian of the transformation 

^ dx'^' 

Assume now, that the metric components are only Lipschitz continuous with re- 
spect to across S. Then differentiating (jS.ip in the direction w = w^-^s^ we obtain 

[gaf^M^^ = Vpig^uA^" + g^.uJUiJU^" + 9^.uW!i,.]w'' , (5.2) 

where [/] denotes the jump in the quantity / across the shock surface S. Thus, since 
both g and Ja are in general Lipschitz continuous across S, the jumps appear only on 
the derivatives. Equation (|5.2p gives a necessary and sufficient condition for the metric 
g to be C^'^ in x" coordinates. Namely, taking w = (|5.2p implies that [50/3,7] = 
for every a, 7 = 0, 3 if and only if 

[J&Ma^.. + [Jl^Vlig^.. + Jlir^[9f.u,j] = . (5.3) 

Note that if the coordinate transformation is C^, so that Ja is C^, then the jumps 
in J vanish, and ()5.2p reduces to 

which is tensorial because the non-tensorial terms cancel out in the jump [(7^/3, -y]. Since 
tensor transformations preserve the zero tensor, it is precisely the lack of covariance 
in (j5.2p for C^'^ transformations that provides the necessary degrees of freedom, (the 
jumps [Ja,'y] in the first derivatives of the Jacobian), that make it possible for a Lips- 
chitz metric to be smoothed by coordinate transformation at points on a single shock 
surface, illustrating that there is no hope of lifting the metric regularity by coordinate 
transformations that are C^. 

Equation (|5.3p is linear in the jumps in the derivatives of the Jacobians, and our 
intention is to use this to solve for the [Ja,-y] associated with a given C^'^ coordinate 
transformation. To this end, suppose we are given a single radial shock surface S in 
SSC locally parameterized by 

j:{t,e,(i)) = {t,x{t),e,^). (5.4) 
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For such a hypersurface in Standard Schwarzschild Coordinates (SSC), the angular 
variables play a passive role, and the essential issue regarding smoothing the metric 
components by C^'^ coordinate transformation, lies in the atlas of {t, r)-coordinate 
transformations. Thus we restrict to the atlas of (i, r)-coordinate transformations for 
a general C^'^ metric in SSC, c.f. (j2.10p . The following lemma gives the unique solution 
[Ja,-/] of (|5.3p for (t, r)-transformations of C^'^ metrics g in SSC. 

Lemma 5.1. Let 

g^^ = -A(t, r)dt^ + B{t, r)dr^ + r^dn"^ , 

be a given metric expressed in SSC, let T, denote a single radial shock surface (15. 4|) 
across which g is only Lipschitz continuous. Then the unique solution [Ja,'y] of (|5.3|) 
which satisfies the integrability condition, {c.f ()2.3p ). 

= > (5-5) 

is given by: 



Wo,tl — 2 \ A A ' ~ 2 \ A A 



[^1,] = "2 [-B-Ji + ^-^1 J . [-^1,.] = "2 ^-B--^! + -g--^! J • (^-6) 

(H^e use the notation & {tj'''} ^'^^ a,/5 € {0,1}, so i/iai t,r are used to denote 
indices whenever they appear on the Jacobian J.) 

Proof: Equation ()5.3p as an inhomogeneous 6x6 linear system in eight unknowns 
[Ja,'y]- Imposing the integrability condition in the form of (15. 5p gives two additional 
equations which complete ()5.3p to an 8 x 8 system which is uniquely solvable for [^0,7]. 
The result is a purely algebraic system whose unique solution (j5.6[) we obtain by a 
lengthy calculation aided by MAPLE, (c.f. [8] for details.) □ 

Condition (|5.3p is a necessary and sufficient condition for [^a/j,^] = at a point on 
a smooth single shock surface. Because Lemma |5. II tells us that we can uniquely solve 
()5.3p for the Jacobian derivatives, it follows that a necessary and sufficient condition for 
[50^3,7] = is also that the jumps in the Jacobian derivatives be exactly the functions 
of the jumps in the original SSC metric components recorded in ()5.6p . In light of this. 
Lemma [5T] immediately implies the following corollary: 

Corollary 5.2. Let p be a point on a single smooth shock curve 7, and let g^^^ be a 
metric tensor in SSC, which is C^'"^ across 7 in the sense of Definition Suppose 
Ja is the Jacobian of an actual coordinate transformation defined on a neighborhood 
Af of p. Then the metric in the new coordinates g^p is in C^'^{M) if and only if Ja 
satisfies (iSlBllFI 



-'^'-'Note that to lift the metric regularity, the Jacobian must mirror the regularity of the metric in 
order to compensate for all discontinuous first order derivatives of the metric by its own discontinu- 
ous first order derivatives. This explains why only C^'^ transformations can possibly lift the metric 
regularity from C"'^ to C^'^ , and C^'" does not suffice for a ^ 1. 
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We conclude that (j5.6p is a necessary and sufficient condition for a coordinate trans- 
formation to lift the regularity of an SSC metric from C^'^ to C^'^ at a point on a single 
smooth shock surface. The condition relates the jumps in the derivatives of the Ja- 
cobian to the jumps in the metric derivatives across the shock. This establishes the 
rather remarkable result that there is no algebraic obstruction to lifting the regularity 
in the sense that the jumps in the Jacobian derivatives can be uniquely solved for in 
terms of the jumps in the metric derivatives, precisely when the integrability condition 
()5.5p is imposed. The condition is a statement purely about spherically symmetric 
spacetime metrics in SSC coordinates because neither the RH conditions nor the Ein- 
stein equations have yet been imposed. But we know by Israel's theorem that the RH 
conditions must be imposed to conclude that smoothing transformations exist. The 
point then, is that to prove the existence of coordinate transformations that lift the 
regularity of SSC metrics to C^'^ at p G S, we must prove that there exists a set of 
functions Ja defined in a neighborhood of p, such that (j5.6p holds at p, and such that 
the integrability condition ()2.3p . (required for Ja to be the Jacobian of a coordinate 
transformation), holds in a whole neighborhood containing p. In the next section we 
give an alternative proof of Israel's Theorem by showing that such Ja always exist in a 
neighborhood of a point p on a smooth single shock surface, and the following section 
we prove that no such functions exist in a neighborhood of a point p of shock wave 
interaction, unless Det {Ja) = at p. 

6. Metric Smoothing on Single Shock Surfaces and a Constructive 

Proof of Israel's Theorem 

We have shown in Corollary 15.21 that (j5.6p is a necessary and sufficient condition on a 
Jacobian derivative Ja for lifting the SSC metric regularity to C^'^ in a neighborhood 
of a shock curve. We now address the issue of how to obtain such Jacobians of actual 
coordinate transformations defined on a whole neighborhood of a shock surface. For 
this we need to find a set of functions Ja that satisfies (|5.6p . and also satisfies the 
integrability condition (j2.3p in a whole neighborhood. In this section we show that this 
can be accomplished in the case of single shock surfaces, thereby giving an alternative 
constructive proof of Israel's Theorem for spherically symmetric spacetimes: 

Theorem 6.1. (Israel's Theorem) Suppose g^y is an SSC metric that is C^'^ across 
a radial shock surface 7 in the sense of Definition \4-i\ such that it solves the Einstein 
equations ()2.1ip - (j2.14p strongly away from 7, and assume T^'^ is everywhere bounded 
and in away from 7. Then around each point p on j there exists a C^'^ coordinate 
transformation of the {t,r)-plane, defined in a neighborhood J\f of p, such that the 
transformed metric components gap are C"^'^ functions of the new coordinates, if and 
only if the RH jump conditions ()3.4p . ()3.5p hold on ^ in a neighborhood of p. 

The main step is to construct Jacobians acting on the (t, r)-plane that satisfy the 
smoothing condition (j5.6p on the shock curve, the condition that guarantees [^0/3,7] = 0. 
The following lemma gives an explicit formula for functions Ja satisfying (15. 6p . The 
main point is that, in the case of single shock curves, both the RH jump conditions 
and the Einstein equations are necessary and sufficient for such functions J a to exist. 

Lemma 6.2. Let p he a point on a single shock curve 7 across which the SSC metric 
g^u is Lipschitz continuous in the sense of Definition \4-l\ in a neighborhood M of p. 
Then there exists a set of functions Ja € C^'^{M) satisfying the smoothing condition 
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()5.6|) on 'y n M if and only if (j4.8p holds on 7 n M. Furthermore, all Ja that are in 
C^'^{J\f) and satisfy (j5.6p on j r\J\f are given by 



4{t,r) 
J*(t,r) 
4it,r) 



[Ar]cPit) + [BMt) 
4A o 7(t) 

+ mat) 

4Ao^{t) 
[BtW) + [^,]a;(t) 
4Bo7(t) 





- r + <I)(t,r) 




-r| + A^(t,r) 




- r| + n{t,r) 


x(t) 


— r + r) 



/or arbitrary functions Z, N G C'^'"'^(A/'), where 

= <|)o7, a; = r2o7, 1/ = N o Q = Z o ^ . (6-2) 
Moreover, each arbitrary function U = ^,0,, Z or N satisfies 

[Ur] = = [Ut]. (6.3) 



Proof: Suppose there exists a set of functions Jq € C^'^{J\f) satisfying (j5.6p . then 
their continuity imphes that tangential derivatives along 7 match across 7, that is 

[<J = -i[J^,r] (6.4) 
for all ^ G and q € {0, 1}. Imposing (I6.4p in (15.61) and using (14. 3|) - (14. 4|) yields 

(1^:81). 

To prove the opposite direction it suffices to show that all t and r derivatives of Ja, 
defined in (j6.ip . satisfy (15. 6p for all /i € and a € {0,1}. This follows directly 

from (14. Sp . (|4.4p and (j4.8p . upon noting that (j6.2p implies the identities 

= J* o 7, 1/ = J* o 7, (J = o 7, ( = J[ o 7 . (6.5) 

This proves the existence of functions Jq satisfying (|5.6p . Applying (the one shock 
version of) Corollary 14.31 (which allows $ to have the lower regularity $ G C^'^ but 
imposes the jumps (j7.4p along 7), confirms that all such functions can be written in 
the canonical form (|6.ip . □ 

To complete the proof of Israel's Theorem, we must prove the existence of coordinate 
transformations — >■ that lift the C^'^ regularity of g^i,^ to C^'^. It remains, then, 
to show that the functions Ja defined above in ansatz (j6.ip can be integrated to 
coordinate functions, i.e., that they satisfy the integrability condition ()2.3p in a whole 
neighborhood. This is accomplished in the following two lemmas. 

Lemma 6.3. The functions Ja defined in (j6.ip satisfy the integrability condition \2. 3j) 
if an only if the free functions <1>, fi, N and Z satisfy the following system of two PDE's: 

ia\X\ + <^t) (/3|X1 + N) + <i>r ie\X\ + Z) - {a\X\ + $) (/3|X| + Nt^ (6.6) 

-Nr{6\X\+n) + fH{X) = 

(^6\X\ + Qt^ il3\X\ + N) + Qr {e\X\ + Z) - {e\X\ + Zt) {a\X\ + $) (6.7) 

-Zri5\X\+n) + hH{X) = 0, 
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where X{t,r) = x{t) — r, H{-) denotes the Heaviside step function, 

[Ar]m + [BtMt) . [AMt)±mm. 

" 4Ao-f{t) ' ^ 4Ao'y{t) 

^ [Bt]<P{t) + [Br]uj{t) _ _ [BMt) + [Br]C{t) _ 

ABoj{t) ' ^ ABoj{t) 

and 



(6.8) 



/ = {/36-a€)\X\ + axN-l3x<^> + f3n-aZ 
h = {^6 -ae)x\X\ + 6xN -ex^ + en-6Z , (6.9) 
where, a, 6 and e are functions of t and f and h are in 

The proof of Lemma 16.31 follows by substituting ansatz (|6.ip into the integrability 
condition ()2.3p and identifying the terms in the resulting first order differential equa- 
tions for Ja. (For details see [8].) 

The proof of Israel's Theorem is complete once we prove the existence of solutions 
^,n,N and Z of ([62]), dSZl) that are C°'\ such that they satisfy dOD- For this 
it suffices to choose N and Z arbitrarily, so that ()6.6p . ()6.7p reduces to a system of 
2 linear first order PDE's for the unknown functions ^ and i}. The condition ()6.3p 
essentially imposes that fl, N and Z be across the shock 7. Since ()6.6p . ()6.7p are 
linear equations for ^> and Q, they can be solved along characteristics, and so the only 
obstacle to solutions <^ and Q with the requisite smoothness to satisfy the condition 
(16. 3p . is the presence of the Heaviside function H{X) on the right hand side of (16. 6p . 
(16. 7p . Lemma [Ql thus isolates the discontinuous behavior of equations (16. 6p . ()6.7p in 
the functions / and /i, the coefficients of H. Israel's theorem is now a consequence 
of the following lemma which states that these coefficients of H{X) vanish precisely 
when the RH jump conditions hold on 7. (See [8] for details.) 

Lemma 6.4. Assume the SSC metric g^y is C^'^ across 7 and solves the first three 
Einstein equations strongly away from 7. Then the coefficients f and g of H{X) in 
(6. 6\) . (6.1) vanish on 7 if and only if the RH jump conditions \2. hold on 7. 



We can now complete the proof of Israel's Theorem. Assuming that the Einstein 
equations hold strongly away from the shock curve (in fact, it suffices to assume that 
only the first three equations hold), we have that there exist functions Ja satisfying 
the smoothing condition (15. 6p if and only if the RH jump conditions hold (c.f. Lemma 
16. 2p . Furthermore, by lemmas 16.31 and 16.41 a solution to the integrability condition 
with the required regularity holds if and only if the RH jump conditions hold (in the 
sense of (14. 8p ). Thus, under the assumption that the Einstein equations hold strongly 
away from 7, we can integrate the Jacobians J a to coordinate functions that smooth 
the metric g to C^'^ if and only if the RH jump conditions hold. This completes the 
proof of Theorem [6l □ 

7. Shock Wave Interactions as Regularity Singularities in GR - 
Transformations in the (t,r)-PLANE 

The main step in the proof of Theorem 11.11 is to prove that there do not exist C^'^ 
coordinate transformations of the (t, r)-plane in a neighborhood of a point p of regular 
shock wave interaction in SSC that lifts the regularity of the metric g from C^'^ to C^'^ 
in a neighborhood of p. We then prove in Section [8] that no such transformation can 
exist within the full C^'^ atlas that transforms all four variables of the spacetime, i.e.. 
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including the angular variables. We formulate the main step precisely for lower shock 
wave interactions in in the following theorem, which is the topic of this section. 
A corresponding result applies to upper shock wave interactions in i?^, as well as two 
wave interactions in a whole neighborhood of p. 

Theorem 7.1. Suppose that p is a point of regular shock wave interaction in SSC, 
in the sense of Definition 13. il for the SSC metric g^y. Then there does not exist a 
C^'^ coordinate transformation x"o(x^)^-'^ of the {t,r)-plane, defined on A/'flM^ for a 
neighborhood M of p in R^, such that the metric components gap are functions of 
the coordinates in AfnM."^ and such that the metric has a non-vanishing determinant 
at p, (that is, such that lim Det (gaB(q)) 7^ 0)0 

In the remainder of this section we outline the proof of Theorem 17. H which mir- 
rors the constructive proof of Israel's Theorem 16.11 in that it uses the extension ()7.ip 
of ansatz ()6.ip to construct all C^'^ coordinate transformations that can smooth the 
gravitational metric to C^'^ in a neighborhood of a point p of regular shock wave inter- 
action. The negative conclusion is then reached by proving that any such coordinate 
transformation must have a vanishing Jacobian determinant at p. But now, to prove 
non-existence, we must show the ansatz ()7.ip is general enough to include all C"'^ 
Jacobians that could possibly lift the regularity of the metric. For this we use condi- 
tion ()5.6p to construct a canonical form for the Jacobians in a neighborhood of p, that 
generalizes ()6.ip to the case of two shock curves, with the weaker assumption of C^'^ 
regularity on the functions $, ri, Z, N. We conclude the proof by showing that this 
canonical form is inconsistent with the assumption that Det [gap) 7^ at p, by using 
the continuity of the Jacobians up to p. 

To implement these ideas, the main step is to show that the canonical form 14.121 
of Corollary 14.31 can be applied to the Jacobians Ja in the presence of a shock wave 
interaction. The result is recorded in the following lemma: 

Lemma 7.2. Let p he a point of regular shock wave interaction in SSC in the sense 
of Definition \3. ll corresponding to the SSC metric g^y defined on A/" PI M^. Then there 

exists a set of functions Ja G C°'^(A/' fl M^) satisfying the smoothing condition (15. 6p 
on 7i n AA, i = 1, 2, if and only if (14. Sp holds on each shock curve 7^ n AA. In this case, 
all Ja in C*^'^(A/' n M?_) assume the canonical form 



^ai(t) \xi{t) 



r\ + ^{t,r) 



i 



J{{t,r) 



r| + N{t,r) 



i 



Y^6i{t)\xi{t) 



r\ + Q{t,r) 




(7.1) 



where 



ai{t) 



[Ar\i ^i{t) + [Bt]i UJijt) 

4Ao7i(t) 



Note that Theorem 11.11 states the non-existence of coordinates on an entire neighborhood A/" of 
p in M^, but here we have prove the stronger result that such coordinates do not exist on the upper 
or lower half planes separately. 
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yi{t) + [Bt]i Ci{t) 






[Bt]i 


^^{t) + [Br]^ UJi{t) 




4So7,(t) 


{Bt]i 


Ui{t) + [B,]i Ci{t) 


^Boji{t) 



m 
m -■ 

ei{t) = TJl'^T^'''^' ' (7-2) 

with 

= $o7j, uJi = no'yi, Ci = Zo'yi, Vi = N o -fi , (7.3) 
and where ^,il,Z,N E C^'^{M riB?^) have matching derivatives on each shock curve 

[UAi = = [Ut]i , (7.4) 

forU = <^,n,Z,N, t G (-e,0). 



The essence of the canonical form (j7.ip is that the jumps in derivatives across 
the shock waves have been taken out of the functions ^,Q,Z,N in (j7.4p . We now 
have a canonical form for all functions Jq that meet the necessary and sufficient 
condition (j5.6p for [(7q,/3,^] = 0. However, for Jq to be proper Jacobians that can 
be integrated to a coordinate system, we must use the free functions <I>, Q, Z, N to 
meet the integrability condition (j2.3p . To finish the proof of Theorem 17. H we show 
that, as a consequence of (|7.4p . (that is, the free functions are regular at the shocks), 
the Jacobian determinant Det Ja must vanish at the point of shock interaction, which 
then implies Det [gap] = 0. 

Thus, using the canonical form (j7.ip restricted to the shock curve and taking the 
determinant of the resulting Ja leads directly to 

Det (JH o 7,(t)) = ( J* J[ - J* J^) l^^(i) = Mt)Qit) - ly^{t)uJ^{t). (7.5) 

Since Ja is continuous, we obtain the same limit t ^ for i = 1,2, 

lim Det (J^ o 7,(t)) = UO)UO) - z^i(0)wi(0) = 0oCo - '^o^o- (7.6) 

t-i>0+ 

Therefore, the final step in the proof of Theorem 17.11 is the following lemma: 

Lemma 7.3. Let p ^ Af be a point of regular shock wave interaction in SSC in the 
sense of Definition \3.1[ Then if the integrability condition 

holds in A/'nM?_ for the functions Ja defined in JZiD , (so that Q,, N and Z satisfy 
dZaD), then 

Ib + ^) - ^2) (0oCo - z^o^o) = 0. (7.8) 

Proof: Substituting the J^ in ([7l]) into dLT]) gives equations (16. 6p . (16. 7p except that 
we now sum over two shock curves instead of one. The difference is the appearance of 
additional mixed terms in the coefficients / and h of the discontinuous terms multiply- 
ing the Heaviside function H{X). The proof is accomplished by showing that, unlike 
/ and g in (16. 6p . (16. 7p . these mixed terms do not vanish by the jump conditions for 
the Einstein equations alone. Finally, a lengthy calculation to evaluate the limit t ^ 
demonstrates that imposing the condition that these additional mixed terms should 
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vanish, which is necessary for (|7.4p to hold, imphes the final equation (j7.8p . (See ^Sj 
for details.) □ 

To finish the proof of Theorem 17.11 observe that the first three terms in (j7.8p are 
nonzero by our assumption that shock curves are non-null, and have distinct speeds 
at t = 0. Thus (|7.8p implies 

Det J^{p) = ((/>oCo - i^ou^o) = 0, (7.9) 

as claimed. □ 

In summary, we remark that at first there appears to be more than enough freedom 
to choose the free functions ^,0,, Z,N in the canonical form of Lemma [721 to arrange 
for the discontinuous term in the integrability condition to vanish, (just as in Lemma 
16.21 the main step leading to Israel's Theorem). This together with the fact that the 
derivatives of Ja are uniquely solvable in condition (|5.6p . lead us to believe until the 
very end that one could construct coordinates in which ga/s was C^'^. But at the very 
last step, taking the limit to the point p of shock wave interaction, we find that the 
condition (|7.4p . expressing that [(7«/3,7] vanishes at shocks, has the effect of freezing 
out all the freedom in Z, N , thereby forcing condition ()7.9p . implying that the 
determinant of the Jacobian must vanish at p. The answer was not apparent until the 
very last step, and thus we find the result quite remarkable and surprising. 

8. Shock Wave Interactions as Regularity Singularities in GR - 

THE Full Atlas 

For the proof of Theorem ll.il we have established the nonexistence of C^'^ coordi- 
nate transformations in the (t, r)-plane that can map a C^'^ regular SSC metric g^;^ 
over to a C^'^ metric ga/s- It remains to extend this result to the full atlas of coordi- 
nate transformations that depend on all four coordinate variables, including the SSC 
angular variables. In this section we outline the proof, (see [8j for details). 

So assume for contradiction there exist coordinates in which the metric gij is C^. 
In general gij is not of the box diagonal form (|2.9p . however, one can always transform 
back to a metric in box diagonal form (c.f. fi5\ chapter 13). The point is now that 
this transformation is in C^, since Killing's equation 

yields a regular Killing vector field X, provided gij is in C^, which if integrated up 
to coordinates yields a coordinate transformation. Therefore the resulting metric 
in box diagonal form gap is regular. Now ga(3 can be taken over to SSC by a 
C^'^ coordinate transformation, which, together with its inverse, act only on the {t, r) 
variables. This contradicts Theorem 17. 1^ and completes the proof of Theorem 11.11 

9. The Loss of Locally Inertial Frames 

Finally we discuss the non-existence of locally inertial frames around a point of 
regular shock wave interaction. This is in surprising contrast to the case of points 
on single shock surfaces for which locally inertial coordinate frames always exist, (c.f. 
|llj). To start, we clarify what we mean by a locally inertial frame: 

Definition 9.1. Let p be a point in a Lorentz manifold and let he a coordinate 
system defined in a neighborhood of p. We call x^ a locally inertial frame around p if 
the metric gij in those coordinates satisfies: 
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(i) 9ij{p) = Vij, Vij = diag{-l, 1, 1, 1), 
(a) 9ij,l{p) = ^ for alliJJ G {0, ...,3}, 

(in) Qij^kl ore hounded in every compact neighborhood of p. 

When the metric components gij satisfy (l)-(3) at point p, we say gij is locally 
Minkowskian (or locahy flat or locahy inertial) at p. By Theorem II. H there exist 
second order derivatives of the metric which are pointwise unbounded in every neigh- 
borhood of p. Therefore, the fohowing Corohary is a straightforward consequence of 
Theorem 11.11 

Corollary 9.2. Let p be a point of regular shock wave interaction in SSC in the sense 
of Definition \3.1\ Then there does not exist a C^'^ coordinate transformation such that 
the resulting metric gij is locally Minkowskian around p. 

10. Conclusion 

Our results show that points of shock wave interaction give rise to a new kind of 
singularity which is different from the well known singularities of GR. The famous 
examples of singularities in GR are either non-removable singularities beyond phys- 
ical spacetime, (for example the center of the Schwarzschild and Kerr metrics, and 
the Big Bang singularity in cosmology where the curvature cannot be bounded), or 
else they are removable in the sense that they can be transformed to locally iner- 
tial points of a regular spacetime under coordinate transformation, (for example, the 
apparent singularity at the Schwarzschild radius, the interface at vacuum in the in- 
terior Schwarzschild, Oppenheimer-Snyder [9j, Smoller- Temple shock wave solutions 
|13t I12j. and any apparent singularity at smooth shock surfaces that become regular- 
ized by Israel's Theorem, [Hill]). In contrast, points of regular shock wave interaction 
are non-removable singularities that propagate in physically meaningful spacetimes in 
GR, such that the curvature is uniformly bounded, but the spacetime is essentially not 
locally inertial at the singularity. For this reason we call these regularity singularities. 

Since the gravitational metric tensor is not locally inertial at points of shock wave in- 
teraction, it begs the question as to whether there are general relativistic gravitational 
effects at points of shock wave interaction that cannot be predicted from the compress- 
ible Euler equations in special relativity alone. Indeed, even if there are dissipativity 
terms, like those of the Navier Stokes equationsj^ which regularize the gravitational 
metric at points of shock wave interaction, our results assert that the steep gradients 
in the derivative of the metric tensor at small viscosity cannot be removed uniformly 
while keeping the metric determinant uniformly bounded away from zero, so one would 
expect the general relativistic effects at points of shock wave interaction to persist. We 
thus wonder whether shock wave interactions might provide a physical regime where 
new general relativistic effects might be observed. Said differently, a regularity singu- 
larity is not hidden behind an event horizon, so it is a sort of counterexample to the 

"'^■^This condition ensures that the physical equations in GR (which are tensorial) difTer from the 
corresponding equations in flat Minkowski space by only gravitational effects, i.e., effects that are 
second order in the metric derivatives. In most of the literature on GR, (c.f. [4]), the gravitational 
metric is assumed to be at least C^'^, (c.f. [3]), which then directly implies condition (3) of Definition 
19.11 At this stage it is not clear to the authors whether or not there exist coordinates that could play 
a (physically satisfying) role as locally inertial frames by satisfying (l)-(2), but not (3). 

^■^The issue of how to incorporate a relativistic viscosity that meets the speed of light bound is 
problematic, |15) . 



REGULARITY SINGULARITIES IN GENERAL RELATIVITY 



21 



cosmic censorship conjecture in the sense that it gives rise to unbounded second order 
metric derivatives, which by themselves might yield physically measurable effects that 
resemble some effects of unbounded curvature. 
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